Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 15 August 2008 (MN WT^ style file v2.2) 

Time-Dependent of Accretion Flow with Toroidal 
Magnetic Field 



00 

o 
o 

(N 

< 



6 



(N 
> 

en 

00 

a^ 
o 

oo 
O 



Alireza Khesali * and Kazem Faghei f 

Department of Physics, Mazandaran University, Babolsar, Iran 



ABSTRACT 

In the present study time evolution of quasi-spherical polytropic accretion flow with 
toroidal magnetic field is investigated. The study especially focused the astrophysically 
important case in which the adiabatic exponent 7 = 5/3. In this scenario, it was 
assumed that the angular momentum transport is due to viscous turbulence and used 
a-prescription for kinematic coefficient of viscosity. The equations of accretion flow are 
solved in a simplifled one-dimensional model that neglects the latitudinal dependence 
of the flow. In order to solve the integrated equations which govern the dynamical 
behavior of the accretion flow, self-similar solution was used. The solution provides 
some insight into the dynamics of quasi-spherical accretion flow and avoids many of 
the strictures of the steady self-similar solution. The effect of the toroidal magnetic 
field is considered with additional variable f3[— Pmag/Pgas], where Pmag and Pgas are 
the magnetic and gas pressure, respectively. The solution indicates a transonic point in 
the accretion flow, that this point approaches to central object by adding strength of 
the magnetic field. Also, by adding strength of the magnetic field, the radial-thickness 
of the disk decreases and the disk compresses. It was analytically indicated that the 
radial velocity is only a function of Alfv'en velocity. The model implies that the flow 
has differential rotation and is sub-Keplerian at all radii. 
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1 INTRODUCTION 

Accretion is the main source of energy in many astrophysi- 
cal objects including different types of binary stars, binary 
X-ray sources, quasars, and Active Galactic Nuclei (AGN). 
Tliough the first development of accretion tlieory started 
a long time ago (Bondi & Hoyle 1944; Bondi 1952), inten- 
sive development of the theory began after the discovery of 
the first X-ray sources (Giacconi et al. 1962) and quasars 
(Schmidt 1963). Since the removal of the angular momen- 
tum process operates on slower timescales as compared to 
free-fall time, the infalling gas with sufficiently high angular 
momentum can form a disklike structure around a central 
gravitating body that can be thin or thick depending upon 
their geometrical shapes. The models of thin accretion disks 
are perhaps better developed and seem to have good ob- 
servational bcisis (Shakura & Sunyaev 1973). However, for 
the thick accretion disks, no fully developed model exists, 
and there remain many theoretical uncertainties about their 
structure and stability (Banerjee et al. 1995; Ghanbari & 
Abbassi 2004; Ghanbari et al. 2007). 
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It is thought that accretion disks, whether in star- 
forming regions, in X-ray binaries, in cataclysmic variables, 
or in the centers of active galactic nuclei, are likely to be 
threaded by magnetic fields. Consequently, the role of mag- 
netic fields has been analyzed in detail by a number of in- 
vestigators (Blandford & Znajek 1977; Lubow et al. 1994; 
Banerjee et al. 1995; Shadmehri 2004) .A mechanism for an- 
gular momentum transport is another key ingredient in the- 
ory of accretion processes and still many theoretical uncer- 
tainties remain about its nature. As originally pointed out in 
Lynden-Bell (1969) and Shakura & Sunyaev (1973), a mag- 
netic field can also contribute to the angular momentum 
transport. A robust mechanism of the excitation of magne- 
tohydrodynamical (MHD) turbulence was shown to operate 
in accretion disks due to the Magneto- Rotational Instability 
(MRI) (Balbus & Hawley 1998; Machida et al. 1999; Begeb 
man & Pringle 2007). 

The toroidal magnetic fields have been observed in the 
outer regions of YSO discs ((Aitken et al. 1993; Wright et al. 
1993; Greaves et al. 1997) and in the Galactic center (Novak 
et al. 2003; Chuss et al. 2003). Accretion disks containing 
toroidal magnetic field have been studied by several authors 
(Fukue & Okada 1990; Geroyannis & Sidiras 1992,1993,1995; 
Banerjee et al. 1995; Terquem & Papaloizou 1996; Machida 
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et al. 1999; Liffman & Bardou 1999; Rempel 2006; Begel- 
man & Pringle 2007; Akizuki and Fukue 2006). Fukue and 
Okada (1990) examined the oscillations of a gaseous disk 
which were penetrated by toroidal magnetic fields. Geroy- 
annis and Sidiras (1992,1993) described differentially rotat- 
ing polytropic models distorted by toroidal magnetic fields. 
Also, Geroyannis and Sidiras (1995) considered dissipative 
effects by viscous friction of differentially rotating visco- 
polytropic models that were further distorted by a toroidal 
magnetic field. Banerjee et al. (1995) presented a toroidal 
magnetic field that was generated by interaction rotating 
plasma and dipolar magnetic field of central object; they 
showed that toroidal magnetic field has an important effect 
in structure of the disk. Terquem and Papaloizou (1996) 
studied the a linear stability of a differentially rotating disk 
containing a purely toroidal magnetic field. They presented 
disks containing a purely toroidal magnetic field are always 
found to be unstable. Machida et al. (1999) considered three- 
dimensional global magnetohydrodynamical simulation of a 
torus treated by toroidal magnetic fields. Akizuki and Fukue 
(2006, hereafter AF) examined the effect of toroidal mag- 
netic field on a viscous gaseous disk around a central ob- 
ject under an advection dominated stage. Assuming steady 
and axisymmetric flow and using steady self-similar method, 
they presented the nature of the disk was significantly dif- 
ferent from that of the weakly magnetized case. 
In this study, we want to explore how the dynamic of a rotat- 
ing and accreting viscous gas depends on its toroidal mag- 
netic field. By solving MHD equations for accreting gases 
that are self-similar in time, we will answer this question. 
We assume that turbulent viscosity is due to angular mo- 
mentum transport of the fluid and there is efficient radia- 
tion cooling in the flow. This paper is organized as follows. 
In section 2, the general problem of constructing a model 
for quasi-spherical magnetized polytropic accretion flow is 
defined. In section 3, self-similar method for solving the in- 
tegrated equations which govern the dynamical behavior of 
the accreting gas is utilized. The summary of the model is 
presented in section 4. 



2 GENERAL FORMULATION 

We use spherical coordinate (r, d, ip) centered on the accret- 
ing object and make the following standard assumptions: 

(i) The accreting gas is a highly ionized gas with infinitive 
conductivity; 

(ii) The magnetic field has only an azimuthal component; 
(iii) The gravitational force on a fluid element is character- 
ized by the Newtonian potential of a point mass,^ = — * , 
with G representing the gravitational constant and M, 
standing for the mass of the central star; 
(iv) The equations written in spherical coordinates are con- 
sidered in the equatorial plane — ^ and terms with any 6 
and (fi dependence are neglected, hence all quantities will be 
expressed in terms of spherical radius r and time t; 

(v) For the sake of simplicity, the self-gravity and general 
relativistic effects have been neglected; 

(vi) The equation of state for the accreting gas is Pgas ~ 
Kp^ with 7 and K being constant. 



The macroscopic behavior of such system can be analyzed 
by perfect magnetohydrodynamics approximation. As stated 
in the introduction, the study focused on analyzing the role 
of toroidal magnetic fleld and viscosity in an accreting gas. 
Thus, the basic equations are the continuity Equation, 



the equations of motion, 
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the polytropic equation. 
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where i}{= v^jr) is the angular speed and v is the kine- 
matic viscosity coefficient. As was mentioned above, our un- 
derstanding of turbulent viscosity is incomplete, and for this 
reason we adopt an empirical prescription, so we employ the 
usual Q-prescription (Shakura & Sunyaev 1973) for the vis- 
cosity which we write in the following form for the kinematic 
coefficient of viscosity. 
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(Narayan & Yi 1994) where a is constant (Tout 2000; King 
et. al. 2007) and Qk is the Keplerian angular velocity, Ke- 
plerian angular velocity is defined by 
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Note that i^ is a function of position and time, since Q.k 
depends on r, and p varying by r and t. To study the effect 
of viscosity, a is used as a free parameter. 

Before initiating to solve the equations (l)-(5), it is con- 
venient to non-dimensionalize the equations. So, the dimen- 
sionless variables are introduced according to 
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Under these transformations and with the use of equa- 
tions (4), (6), and (7), equations (1) and (5) do not change, 
but equations (2) and (3) become 
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3 SELF-SIMILAR SOLUTIONS 

3.1 Analysis 

To grasp the physics of the accreting viscous gas in a toroidal 
magnetic field, the technique of self-similar analysis proves 
to be useful. Of course, this method is familiar from its wide 
range of applications in the full set of equations of MHD 
in many research fields of astrophysics. In self-similar for- 
mulation, the various physical quantities are expressed as 
dimensionless functions of a similarity variable, so it lends 
itself to a set of partial differential equations, such as those 
mentioned above, to be transformed into a set of ordinary 
differential equations. A similarity solution, although con- 
stituting only a limited part of problem, is often useful in 
understanding the basic behavior of the system. So, in or- 
der to seek similarity solutions for the above equations, a 
similarity variable rj is introduced as 

and it is assumed that each physical quantity is given by the 
following form: 



pir,t)=f'R{rj) 
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the exponents n, ei, £2, £3, and £4 are constant which 
must be determined. By substituting the equations (12)-(16) 
into equations (1), (5), (10) and (11), the following general 
results are obtained: 

-1. ^4 = - ,, ^ ^, ,(17) 
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The above results imply each physical quantity retain a sim- 
ilar spacial shape as the flow evolves, but the radius of the 
flow increases proportionally to t^''^ . Also time-dependent 
density, the pressure and the toroidal magnetic fleld are 
varying by 7, on the other hand, they are decreasing by 
time for 7 > 1. 

Here, let us seek time-dependent self-similar of the mass 
accretion rate 



5/3, time-dependent behavior of this quantity will be ap- 
plied in next section. 

Solving equations (1), (10), (11), and (19) under trans- 
formations (12)-(15) in nonmagnetically state, makes it clear 
that behavior of physical quantities in the nonmagnetically 
and the magnetically disk are the same. The result is one of 
the strictures of time-dependent self-similar solution. 

Subsequently, the equations (1), (5), (10), and (11) for 
the dependence of the physical quantities on the similarity 
variables are written as: 
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This is a system of non-linear ordinary differential equations. 
Once a and 7 are selected, the set of equations (22)-(25) can 
be solved. Before solving above equations numerically, it was 
found out that equations (22) and (25) imply 
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where C is constant of integration, and will be calculated in 
the next section. We can rewrite equation (26) in terms of 
the Alfv'en velocity. The Alfv'en velocity in a purely toroidal 
magnetic field is v\ — B'^/Aivp. By using transformations 
of (8), (9), (13), (16), (17), and Alfv'en velocity equation, 
equation (26) can be rewritten in the following form 
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where A = va/v- The result imply that the radial velocity 
of a quasi-spherical accretion flow in presence of toroidal 
magnetic field is a function of Alfv'en velocity. 



M = 



-47rr pVr- 



(19) 



We can non-dimensionalize the equation (19) under trans- 
formation (8) and 



M ^ MM 
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Under transformations (8) and (20), equation (19) does not 
change and its behavior under similarity quantities that are 
implied in equations (12)- (16) can be considered. The simi- 
larity solution shows that the mass accretion rate M is pro- 
portional to t(T^5/-3'/(T^i'. When 7 — 5/3, the mass accre- 
tion rate is independent of time and decreases in 1 < 7 < 



3.2 Inner limit 

When 7 = 5/3, an appropriate asymptotic solution as 77 -^ 
is of the form 
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In order to derive the above relations, the mass accre- 
tion rate and /3 parameter were used, that is ratio of the 
magnetic pressure to the gas pressure. When rj ^ and 
7 = 5/3 the mass accretion rate becomes M ~ — 47ri?oVb, 
and the ratio of the magnetic pressure to the gas pressure 
becomes /3(r;) ~ PoV > where /3o = Bq/Rq . These rela- 
tions were applied to derive equations (32)- (35). In order 
to present of importance of magnetic field in the disk, /3o 
parameter will be used. 

Asymptotic solution shows that a parameter is effective 
in inner edge of the disk and physical quantities are sensi- 
tive to it, i.e., the radial infall velocity increases by adding 
a, the angular velocity is sub-Keplerian for all values of a, 
and the density and the toroidal magnetic field in the inner 
edge of the disk decrease with increasing a. These results 
that are achieved for inner edge of the disk are qualitatively 
consistent with the results of AF. Now, it is possible to de- 
rive approximate constant of integration C in equation (26), 
by using equations (28)- (35) 
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3.3 Numerical solution 

If the value of rii„ is guessed, that is a point very near of the 
center, the equations can be integrated from this point to 
outward through the use of the above expansions. Examples 
of such solutions are presented in Figs 1 and 2. The profiles 
in Fig 1 and 2 are plotted for different l3in that is amount of 
f3 in rjin. From guessed rjin and /3i„, and equations (28), (31), 
(32), and (35), we can achieve /3o = Pin/rj^^ . The delineated 
quantities {rf^'^R, ri^'^V , ...) in Figs 1 and 2 are constant at 
steady self-similar solutions (Narayan & Yi 1994; Narayan & 
Yi 1995; Shadmehri 2004; Akizuki & Fukue 2006; Ghanbari 
et. al. 2007). 

By increasing the /3 parameter, which indicates the role 
of magnetic field in the dynamics of accretion disks, the 
radial-thickness of the disk decreases; the equations (13) 
and (17) imply that compression increases by time. Liff- 
man & Bardou (1999) and Campbell & Heptinstall (1998) 
showed compression of disk in height direction by effect of 
toroidal magnetic field, but they did not consider the effect 
of toroidal magnetic field in the radial-thickness of the disk. 
Also, by adding the /3 parameter, the radial infall velocity 
increases; such property is qualitatively consistent with AF. 
This is due to the magnetic tension terms, which dominate 
the magnetic pressure term in the radial momentum equa- 
tion that assist the radial infall motion. The flow is differen- 
tially rotating, although it is highly sub-Keplerian at large 
radii. 

For investigating of existence of transonic point, the 
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Figure 1. Time-dependent self-similar solution for 7 = 5/3, o = 
0.5, and M = 1.0. The solid lines represent /3i„ = 0.5, the dotted 
lines represent /3i„ = 1.0, and the dashed lines represent /3i„ = 1.5 
that f3i„ is value of f3 in rii„ . 
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Figure 2. Timo-dcpendent self-similar solution for 7 = 5/3, a = 
0.5, and M = 1.0. The solid lines represent /3i„ = 0.5, the dotted 
lines represent f3i„ = 1.0, and the dashed lines represent I3i„ = 1.5 
that /3i„ is value of /3 in rii„. 



square of the sound velocity is introduced that subsequently 
can be expressed as 
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Here, 5* = (jR'^^'^] ' the adiabatic sound velocity in self- 
similar flow, which is rescaled in the course of time. The 
Mach number referred to the reference frame is defined as 
(Fukue 1984; Gaffet & Fukue 1983) 
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is the velocity of the reference frame which is moving out- 
ward here as time goes by. The Mach number introduced so 
far, represents the instantaneous and local Mach number of 
the unsteady self-similar flow. As we see in Fig 2, there is 
a transonic point, that denotes the square of Mach number 
is equal to unit {y? — 1). By adding strength of the mag- 
netic field, transonic point approaches to central object. The 
solution shows that /3 parameter varies by radii and is im- 
portant at larger radii, while in steady self similar solution 
this parameter is constant. The /3 parameter shows that the 
dominate pressure in the outer region of disk is magnetic 
pressure, that this result is consistent with observed YSO 



disks (Aitken et al. 1993; Wright et al. 1993; Greaves, Hol- 
land & Ward-Thompson 1997). 



4 SUMMARY 

In this paper, the equations of time- dependent quasi- 
spherical accretion fiow with toroidal magnetic field have 
been solved by semi-analytical similarity methods. The flow 
is able to radiate efflciency, so we substituted the polytropic 
equation instead energy equation. A solution was found for 
the important case 7 = 5/3 that has differential rotation 
and viscous dissipation. The flow avoids many of the stric- 
tures of steady self-similar solutions (Narayan & Yi 1994; 
Akizuki & Fukue 2006) . Thus, the radial-dependence of cal- 
culated physical quantities in this sense are different from 
steady self-similar solution. 

The fiow has differential rotation in small radii and 
has Keplerian behavior at large radii, that at large radii 
is similar to steady self-similar solutions. Also, The fiow is 
sub-Keplerian at all radii that is consistent with AF when 
they considered disk in moderate strength of the magnetic 
field. The solution shows that in time-dependent of quasi- 
spherical accretion fiow, there is a transonic point, where the 
point approaches to central object by increasing strength of 
the toroidal magnetic field. By increasing strength of the 
toroidal magnetic field, the radial thickness of the disk de- 
creases and disk becomes compress. 

Here, latitudinal dependence of physical quantities is 
ignored, although some authors showed that latitudinal de- 
pendence is important in structure of a disk (Narayan & 
Yi 1995; Ghanbari et. al. 2007). One can investigate lat- 
itudinal behavior of such disks. Also, it is assumed that 
there is efficient radiation cooling in the flow and used poly- 
tropic equation for energy equation. During recent years one 
type of accretion disks has been studied, in which the en- 
ergy released through viscous processes in the disk may be 
trapped within the accreting gas. This kind of flow is known 
as advection-dominated accretion fiow (ADAF). Solution of 
AF shows that physical quantities of the disk vary by advec- 
tion parameter. In future studies, we are going to improve 
our model with a realistic energy equation. 
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